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ABSTRACT 

We propose a new mechanism which explains the existence of enormously sharp edges 
in the rings of Saturn. This mechanism is based on the synchronization phenomenon 
due to which the epicycle rotational phases of particles in the ring, under certain 
conditions, become synchronized with the phase of external satellite, e.g. with the 
phase of Mimas in the case of the outer B ring edge. This synchronization eliminates 
collisions between particles and suppress the diffusion induced by collisions by orders 
of magnitude. The minimum of the diffusion is reached at the center of the synchro- 
nization regime corresponding to the ratio 2:1 between the orbital frequency at the 
edge of B ring and the orbital frequency of Mimas. The synchronization theory gives 
the sharpness of the edge in few tens of meters that is in agreement with available 
observations. 
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1 INTRODUCTION 

Together with very small thickness, extreme sharp edges in 
rings of Saturn are one of t he most outstanding f eature s 



of planetary rings (s e e e.g. FYidman fc Gorkavvil l|l999l ) 



lEspositd (|2002l . |2006| ): iBorderis et all (| 19821 )). Indeed, for 



example the outer B ring edge has a density drop by an or- 
der of magnitude on a distance ~ 10m that is enormously 
sharp compared to the edge distance to Saturn (117580fcm), 
the B ring width (25580fem) and the width of Cassini divi- 
sion (4620fcm). This is especially surprising since the life 
time of the ri ngs is enormously large bein g of about 10^^ or- 
bital periods (|Fridman fc Gorkavvil [l999l ) and the particles 
inside the B ring are quite dense (e.g. there are particles of 
size 10m down to 1cm and smaller wi th a distance between 
them of a bout a few meters and less) (Fridman fc Gorkavvil 
ll999l : lEsposita,2002 , 2006,; Spahn fc Schmid t 2006). Due to 
this about 10 — 100 collisions between particles per orbit 
should wash out all sharp density contrasts in only a few 
orbital periods. 

In many cases, and perhaps always, these sharp edges 
are associated with the gravitational perturbation by a moon 
inside or outside the rings. The abruptness of the tran- 
sition from nearly opaque to practically transparent re- 
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gions was constrai ned by the Voyager Photopolarimeter data 
(iLane et al. I Il982l ) to be smaller than about 100 meters. 
Cassini UVIS occupation profiles show edges sharper than 
tens of meters, and in fact the assumption of a knife-edge 
sharpness is used to constrain the local ring thickness near 
the edge to be on the order of 5 meters only ( (|Colwell et al.l 
|2008| . I2OO9I ) ; ring thickness and sharpness of the edge should 
be of the same order). 



After the pioneering work of IWisdom fc Tremaind 

(| 19881 ) extensive numerical simulations of particle dynam- 
ics have_been_^erformed__by differ en t groups ( see e^ 

Saiol (I1995I): iLewis fc Steward ||2000D: ISeiss et al.l l|2005 



Charnoz et al.l ~ 20071 ) : ISremcevic et all (|2007l )) that allowed 



to establish a number of interesting properties of ring dy- 
namics. However, the problem of sharp ring edges still re- 
mains a mystery. Its solution requires extensive large scale 
numerical simulations with particles of different scales, it 
also may require to go beyond local box simulations invented 
by Wisdom fc Trcmainc (1988). In view of these difficulties 
it seems to be useful to explore certain simplified models that 
bring to surface new qualitative physical effects which can 
be analyzed more directly due to model simplicity. Here, we 
introduce such a simplified model of ring dynamics for outer 
B ring edge called the SYNC model in the following. Numer- 
ical investigations of this model show a striking phenomenon 
of synchronization of the epicycle motions of particles in the 
ring induced, under certain conditions, by a periodical grav- 
itational force of Mimas. In a general context, the synchro- 
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nization phenomenon, which has abundant manifestations in 
scien c e, nature, engi neering and social hfe (|Pikovskv et ahl 
I2OOII : IStrogatdliooa ). can be roughly described as an ad- 
justment of frequencies and phases of oscillators due to in- 
teraction and /or forcing. In the present context, the phases 
of the epicycle rotation of various particles become synchro- 
nized with the phase of periodic gravitational force of Mi- 
mas, and as a result the collisions between particles become 
suppressed by orders of magnitude so that the diffusion in 
the ring also becomes suppressed by orders of magnitude. 
This leads to a maintaince of sharp edges in a certain fre- 
quency range of ratios of epicycle frequency f2 to the Mi- 
mas frequency to. We note that the data of the Cassini mis- 
sion show that the outer B ring edge has the frequency Qb 
which is very close to 2 : 1 resonance with Mimas frequency 
OJ, actually (f^s - ils)/2uj = Q.b/2uj - 1 10"^ - 10"" 
that corresponds to the accuracy in coordinate position of 
about 1 to 10km. Thus f2s is located directly in the mid- 
dle of the synchronization Arnold tongue where the syn- 
chronization effects are the strongest. As a noticeable re- 
mark we mention, that it was Christiaan Huygens who 
discov ered both the Saturn's rin gs (see, e.g. references in 
book (jPridman fc Gorkavvil[l999h ) and the synchronization 
phen omenon (see, e.g. references in book (|Pikovskv et al.l 

The paper is composed as follows: the SYNC model of 
dynamics in the ring is described in Section 2 (with addi- 
tional details given in Appendix A, B), the numerical and 
analytical results are presented in Section 3, the discussion 
of the results is given in Section 4. 



2 DESCRIPTION OF THE SYNC MODEL OF 
RING DYNAMICS 

The SYNC model of the ring dynamics is based on the four 
main ingredients: 

(i) the individual particle dynamics is given by the Hill 
eq uations and is conside r ed in a local box as it is proposed 
by IWisdom fc Tremaind (|l988l ) , the particles are assumed 
to be identical, the dynamics is considered only in the two- 
dimensional plane of the ring; 

(ii) the gravitational force of Mimas is considered as a 
sequence of periodic kicks of fixed amplitude; 

(iii) the collisions between particles are treated on the ba- 
sis of the mesoscop ic multi-particle collision mo del proposed 
by Kapral (see e.g lMalevanets fc Kaprall (120041 ')'): 

(iv) the total energy balance (the process where 
the injection of energy provided by the shear flow 
(|Wisdom fc Tremaind Il988l ) and the gravitational force of 
Mimas is equilibrated by dissipation) is ensured via the 
Nose-Hoover thermostat which is broadly used in molecu- 
lar dynamics sim ulations of lar g e ensembles of interactin g 
particles (see e.g Hooveij (|l999l ') ; iRateitschak et all l|2000l ) ; 
iHoover et all (|2004| )V 



Let us now describe the elements of the model in more 
details: 

(i) The Hill e quations of motion inside a local box 
I'Wisdom fc Trema inc 1988) are 



where Q — \J GMsatum/ is the Kepler frequency of a 
particle of mass nip, Vs is the Kepler shear velocity, Fx(t) 
is the gravitational force of Mimas along axis x directed to 
Saturn. Here Vx , Vy are velocities of local motion in the pres- 
ence of the shear flow. 'With GMsatum = 3.79 x IQ^^m!^ /s^ 
and the radius ao = 1.17 x 10*m we have at the edge of B 
ring Q. = Q.s = 1.52 x 10~*s~^. We normalize all velocities 
by a typical value of epicycle velocity Vep = 0.005m/s that 
gives us equations of motion in a dimensionless form. After 
that the distance is measured in units of a typical epicy- 
cle radius = v^p/Qs ~ 32.7m and time t is replaced by 
Qst- The loca l box has the periodic bound ary conditions as 
those used bv I Wisdom fc Tremaind (|l988h . The size of the 
box is usually taken as a square 5 = 5rs x Sr^. In pres- 
ence of the Nose-Hoover thermostat (see (iv)) and the shear 
flow we found convenient to use the Hamiltonian form of the 
Hill equations as it is described by [stewart (1991;). In this 
formulation the Hamiltonian of the epicycle motion has the 
form Hep = QI = {v^ + Vy)/2 where I is the action of the 
oscillator motion. More details are given in Appendix A. 

(ii) The computation of the field strength Fx is de- 
scribed in Appendix B. Because Mimas passes the local box 
rather fast, the force in dimensionless units has a form of pe- 
riodic delta-function f{t) = Fx{t)/ {rripVep) = tJ2i S{t — lTM) 
with the period tm defined by the dimensionless orbital 
period of Mimas tm = 27rf2s/a; and dimensionless kick 
strength e — 0.64 corresponding to the fixed choice of the 
typical epicycle velocity Vep = 0.005m/s. The force in y- 
direction is neglected since it is much smaller compared to 
the force in a;-direction and gives a small change of Vy com- 
pared to the shear velocity. 

(iii) The collisions of par ticles are performed accord - 
ing to the Kapral algorithm l|Malevanets fc Kaprall [20041 ') . 
Namely, the whole local box S is divided on Ncei cells. Usu- 
ally we use about 100 x 100 cells with the total number of 
particles A'^ = 1000 corresponding to Nep = Trr^N/S ~ 125 
particles inside one epicycle circle and the particle density in 
a cell being 0.1. After a time tk the relative velocities (with 
respect to the motion of the center of mass of the cell) of all 
particles in a given cell are rotated by a random angle. In 
this way the total momentum and energy inside a given cell 
are preserved while the directions of the velocities become 
mixed. It is important to note that during the collision the 
velocities inside the cell are taken as the total physical ve- 
locities of particles, namely x, y including the shear velocity. 
Thus, due to a finite size of the Kapral cells, the shear ve- 
locity always generates appearance of a spreading of local 



velocities Vx 



even if Vx and Vy of two particles coincide 



X = Vx; y = Vy + Vs; Vx = 2Qvy + Fx[t)/mp; 
Vy = -^lvx/2; Vs = -iQ.x/2. 



(1) 



before the "collision" , this is not true for the total velocities 
X, y, and random rotation of the latter leads to the appear- 
ance of non-equal Vx,Vy. In a certain sense the finite size of 
Kapral cells physically acts as a finite size of colliding bod- 
ies. Usually we used tk = 0.5/Qs but the variation of this 
parameter did not affect the main results. 

(iv) In presence of collisions the shear flow and the driv- 
ing Mimas force inject additional energy in the system, that 
is dissipated via different mechanisms, including nonelasitic- 
ity of collisions, interactions with dust, etc. In our simpli- 
fled model, to keep the energy balance, we use the Nose- 
Hoover thermostat which is commonly used for molecular 
dynamics simulations of interacting particles, also in pres- 
ence of external flelds (|Hooveilll99 9l: IRateitschak et al ] |2000l : 
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iHoover et al]|2004l : IChepelianskii et alJlioOTl ). In this ther- 
mostat, which mimics a canonical ensemble, an additional 
friction force acts on a particle i according to 



p, = F. - 7P. , q» = P. ,7 = [{p')/(2mpT) - l]/rH^ (2) 



where piqi are the momentum and coordinate of particle i, 
Fi is an effective "friction" force acting on a particle due 
to collisions and external fields, th is the relaxation time in 
the NH thermostat [usually we used Qsth ~ 16 but we also 
checked that the variation of th does not affect the synchro- 
nization phenomenon (see examples below)] and (p^) means 
the average over all A'^ particles, T = mpvlp/2 is a given tem- 
perature of the thermostat. One can see that the "friction" 
changes the sign with 7, and the variable 7 is driven by 
the deviations of the mean kinetic energy from that at the 
given temperature T, in this way the system is kept near 
this temperature as it should be for the canonical ensemble. 

For the Hamiltonian form of the Hill equations, the fric- 
tion acts only on the action variable that gives in dimension- 
less variables li = —7/^,7 = [(I) — 1]/th'^, where (I) means 
the averaging over all A'' particles (see Appendix A). The 
physical origin of the appearance of such an effective fric- 
tion can be attributed to an average friction force acting on 
a relatively large particle as a result of multiple collisions 
with a dust of small size particles. 

In a ring with extended size distribution the smaller 
particles have in eq uilibrium generally larger dispersion ve- 
locities l|Sald [1993 ). Since the collisions are inelastic, the 
system does however not assume a state of energy equipar- 
tition. Typically, the dispersion velocity of the smallest par- 
ticles is by a factor of several larger than the one of the 
largest particles, depending mildly on the width of the size 
distribution and the inelasticity of the particles. 

At opposition the perturbing moon induces an equal 
excess velocity Av^ to all ring particles, regardless of their 
size. This means that compared to coUisional equilibrium the 
large particles have now a higher excess in random kinetic 
energy than the small ones. In this sense the return to equi- 
librium, mediated by dissipative collisions, affords a cooling 
of the large particles relative to the small ones, and thus, 
an effective friction on large particles. This friction vanishes 
in equilibrium like the Nose-Hoover thermostat. Large par- 
ticles determine the dynamical properties of the ring. 

Another dissipative process could be an ongoing ex- 
change of ring-matter between particles of all size groups 
due to a balance of coagula tion and fragment ation (dynamic 
ephemeral bodies, DEB's, jPavis et al.llQSl )) causing an ef- 
fective dissipation, since the composition and destruction of 
agglomerates are irreversible processes. 

It is interesting to note that the epicycle dynamics is 
rathe r similar to motion of charg ed particles in a magnetic 
field ijFridman fc Gorkavvil [l999l l and due to that there is 
a certain analogy with the synchronization of the Larmor 
motion for two-dimensional electron g as in magnetic and 
microwave fields as it was discussed bv IChepelianskii et al.l 
(|2007l ). The main difference to the present problem is that 
for the electron gas there is no shear. 




Figure 1. Density distribution of particles in the ring in the plane 
of local epicycle velocities (—3 < Vx/vep < 3; —3 < Vy/vep < 3) 
obtained by the numerical simulations with A'^ = 1000 particles 
inside the spacial square box S = 5rs X 5rs where Tb is the epicycle 
radius; the particle density is p = Np/S = iO/r^. The rotation 
frequency ratio is Q/Qg = Q/2ui = 1.15 (left panel) and 1 (right 
panel, synchronized regime); the dimensionless force amplitude of 
Mimas is e = 0.64. The number of Kapral cells is Af^e! = 100 X 
100 = 10"*, the Kapral collisions are done after time tk = 0.5/Os; 
the relaxation time of Nose-Hoover dynamics is Tfj = 16/Qs- The 
data are averaged over time interval ^ 4 W^/fls- Density 
is proportional to color (red/gray for maximum, blue/black for 
minimum). 



3 NUMERICAL RESULTS AND THEIR 
INTERPRETATION 

The results of numerical simulations for the particle density 
distribution in the plane of local epicycle velocities {vx,Vy) 
are shown in Fig.[T]for two values of Q./2uj ratio between the 
epicycle frequency of particles and the double frequency 
of Mimas 2lj. For Q,/2uli — 1.15 the distribution of velocities 
is close to the Maxwell distribution of elliptical form appear- 
ing due to shear and ellipticity of motion given by the Hill 
equations (the distribution becomes close to a symmetric 
one in rescaled velocities Vy = 2vy,Vx = Vx)- The distribu- 
tion is drastically changed for i}/2u) — Q,s/2lu = 1 : almost 
all density is concentrated on a spiral in the velocity plane. 
The physical meaning of this phenomenon is the following: 
for the resonant ratio Q,/2uj = 1 the phases of the epicycle 
rotations become synchronized with the phase of periodic 
Mimas kicks given by fit). While out of the resonance (e.g. 
Q./2uj — 1.15) all epicycle phases are random and indepen- 
dent, in the synchronization regime they are all adjusted to 
the phase of Mimas. As a result, Mimas gives a kick in Vx, 
between kicks the particle velocity decreases doing two rota- 
tions over the spiral, then again it is kicked by Mimas in the 
same state as at previous kick, and so on. Thus the phases of 
this rotational motion are equal to the phase of Mimas and 
equal to each other. This means that the particles rotate in 
synchrony, at each moment of time their rotational veloci- 
ties are equal both in the amplitude and in the direction. 
Therefore, the collisions between them effectively disappear. 
There remains only a residual relative velocity, related to 
the shear velocity and a finite size of the Kapral cells (finite 
size of colliding bodies in the ring) . In the absence of shear 
flow the synchronization can be complete for all par ticles as 
it has been discussed bv IChepelianskii et al.l ()2007l ) for the 
problem of two-dimensional electron gas in magnetic and 
microwave fields. 

Since the collisions are significantly reduced for the syn- 



4 D. L. Shepelyansky, A. S. Pikovsky, J. Schmidt and F. Spahn 




Figure 2. Dependence of the rescaled diffusion rate D/Dq shown 
by color on the rescaled frequency Q/2u! (horizontal axis) and 
driving force strength e (vertical axis) for the range 0.85 ^ 
Q/2uj ^ 1.15 and 5C e 0.7. The color is proportional to D/Dq 
with red/gray for maximum value {D/Dq = 1.26) and blue/black 
for minimum (D/Dq = 4 x 10^*), here Do is the diffusion rate at 
e = 0,n = ns (also D/Dq = 1.7 X 10-5 with Dq = rfn^). Data 
are obtained for time t ^ 10'*/Cs, = 1000, N^^i = 200 X 200, 
other parameters are as in Fig.l. 

chronization regime, the difFusion rate D — {x'^) /t is also re- 
duced by orders of magnitude compared to its typical value 
Do = fjfis- This is clearly seen in Fig. [5] which gives the 
dependence of D on f2 and e. The diffusion suppression takes 
place inside the Arnold tongue where the epicycle rotation is 
synchronized with the Mimas phase. According to the data 
of Fig.[2]the synchronization takes place inside the frequency 
range 

\VL/2u - 1| < se (3) 

with the numerical value of th e constant s ~ 0.08. A ccording 
to the synchronization theory (|Pikovskv et al.|[200ll ) the syn- 
chronization region is given by the dimensionless amplitude 
of the driving force which is 2e/(47rmax(wa;)) that gives s ~ 
0.08 since we have miLX.{vx) ~ 2 (see Fig.[l|. Thus the numer- 
ical dependence is in agreement with the analytical synchro- 
nization theory. The variation of fl is related to the position 
of particle inside the ring. For example, Q,/2u — 1.05 cor- 
responds to the distance Aa; ~ 3ao(fi/2a; — l)/2 ~ QOOOfcm 
from the outer B edge in direction to Saturn. It is interest- 
ing to note that this is of the order of the size of Cassini 
division. 

An interesting property of the relation ([2} is its inde- 
pendence of the relaxation time scale th- Physically, this 
means that th only determines the time scale on which the 
synchronization is reached but it does not affect the domain 
of synchronization. This is in agreement with our numeri- 
cal checks which show that at the Mimas value of e = 0.64 
the synchronization window shows less than 10% variation 
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Figure 3. Dependence of the rescaled diffusion rate D/Dq on 
the rescaled frequency Q/2ui. The full curves are for e = 0.64 
with TVce; = 50 X 50, 100 X 100, 200 X 200 from top to bottom; the 
dashed curve is for e = and Ncei = 100 X 100. Data are obtained 
for time t ^ W^/Qs, N = 1000. Logarithms are decimal. 




Figure 4. Dependence of the synchronization parameter S on the 
rescaled frequency Q/2lj. Parameters are the same as in Fig. 3. 
Logarithms are decimal. 

when (Qsth)'^ varies from 0.1 to 10^^. For {0,sth)^ ~ 1 
this window is increased by about 30%. 

The variation of the diffusion rate D{Q,) when the num- 
ber of Kapral cells Ncei is changed is shown in Fig. |3l 
For Ncet = 50 X 50 the collisions happen rather often and 
there are no signs of synchronization. For Ned = 100 x 100 
the synchronization sets in and the diffusion drops in- 
side the synchronization window. A further increase up to 
Ncei ~ 200 X 200 gives much stronger drop of the dif- 
fusion inside the synchronization window while its size is 
only slightly increased. Outside of this window the diffu- 
sion scales as D cx 1/Ncei- This is rather natural since D 
is proportional to the density of particles inside the cell 
so that D ~ r^r ~ riN / {NceiTK) where F is the effec- 
tive collision rate. We note that in the non-synchronized 
regime the diffusion rate per orbital period is rather large 
being 2TiD/{rlns) « 3 (at N^^i = 100 x 100) and 1 (at 
Ncci = 200 X 200). These values approximately correspond 
to the typical conditions for particles inside B ring. 

Another signature of synchronization can be expressed 
via the synchronization parameter 

S = (vi— Vj)^/(A''^Wep/2). Its dependence on frequency 

is shown in Fig. [3] Inside the synchronization regime S drops 
by almost 4 orders of magnitude. This means that due to 
synchronization the relative collision velocities of particles 
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Figure 5. Dependence of density of particles (in arbitrary units) 
on time Qgt (vertical axis) and position in the ring x/rg (hori- 
zontal axis) for e = 0.6, Q/2lu = 1.1 and zero frequency gradient 
g = 0.. There are A^ce; = 1200 X 120 in the whole space box 
S = 50rs X 5rs; tk = 0.5/fls, th = 4.5/f2s. Initially there are 
A'^ = 60 particles in the left box S = brg X brg and this number 
is kept constant during the computations till the finite moment 
of time ngt = 6.28 X 10* when there are 305 particles in total. 



are very small and therefore the diffusion is also very small. 
At the same time the velocity difference remains finite due to 
finite size of the cells (collision bodies) and the shear flow. In 
absence of the shear flow the collisions disappear completely 
(see also (Che pclia nskii ct a l. 2007)). 

For the chosen value Vep = 0.005m/s we have e = 0.64 
(see Appendix B) that according to ([2]) and the data of 
Figs. I2I3I4I give the synchronization border Q,/2u ~ 1.05. 
This corresponds approximately to xs = 9Q00fcm distance 
from the exact resonance 2:1. The observations give this dis- 
tance to be about xb ~ 1 — lOfcm that is significantly smaller 
than the synchronization border. Of course, the value of Vep 
is not known exactly and may be a factor 2 to 10 different 
from the value chosen above. But even if the actual value is 
10 times larger than still xs remains by 2 orders of magni- 
tude larger than the observed B edge position xb- 

We explain this disagreement in the following way. Tak- 
ing the value e = 0.64 we assume that the particles appear 
initially in the non-synchronized part of the ring, let say with 
Q./2uj ~ 1.2. Due to collisions these particles diffuse closer 
and closer to the synchronization border at Q,/2u ~ 1.05. 
Behind this border the diffusion rate drops D{Q) drastically 
due to the synchronization phenomenon described above. 
This sharp drop of D creates a diffusive shock wave which 
continues to propagate slowly inside the synchronization re- 
gion since there the diffusion remains finite due to the shear 
flow (see Figs. 1213] ). The edge size of this diffusive shock wave 
should be of the size of few epicycle radius Vs since as soon 
as the distance between particles becomes larger than Vs 
the collisions between them completely disappear due to a 
frozen nature of epicycle motion (in a close similarities with 
particles in a magnetic field, see e.g. (|Fridman fc Gorkavvil 
[l999.)). In this way the gradient of density slowly moves in- 
side the synchronization domain keeping the sharpness of 
the edge of a few epicycle radius Va. In the synchroniza- 
tion domain the diffusion is minimal at the exact resonance 
n/2cj — 1 since the synchronization effect is the most strong 



Figure 6. Dependence of density of particles (in arbitrary units) 
on time Hgt (vertical axis in logarithmic scale) and position in 
the ring x/rg (horizontal axis) for e = 0.6, the frequency gradient 
in space is g = 0.002 with Q/2uj = 1.05 at a; = and Q/2lu = 1 
at x/rs =25 corresponding to the outer B ring edge. There are 
Ncei = 1200 X 120 in the whole space box S = 50rg x 5rs (only 
half is shown); tk = 0.5/ fls, th = 4.5/r25. Initially there are 
= 60 particles in the left box S = 5rg X 5rs and this number 
is kept constant during the computations till the finite moment 
of time fist = 6.28 X 10® when there are 305 particles in total. 



there (see Figs. I2I3I4|I . Therefore, the edge of synchronized 
particles spends the most of the time at that place with the 
minimum of diffusion. This is actually the place where the 
outer B ring edge is observed now. 

To give more justification to the above picture we per- 
formed extensive numerical simulations of front propagation 
of particles in the SYNC model. With this aim we intro- 
duced a gradient of the frequency with the distance x so 
that fl{x) = Qo — gx where fio ~ 1.1 is some initial value 
and g ~ 0.002 is the frequency gradient per unit of epicy- 
cle radius in dimensionless units. Examples of the numerical 
simulations are shown in Figs. 5,6. In the absence of the 
gradient, in the non-synchronized regime there is a diffu- 
sive spreading along x as it is clearly seen in Fig. [S] The 
typical diffusive profile x oc Vt is clearly seen. The com- 
putations are done at the fixed constant particle density p 
in the left space box 5rs x 5rs of the total longitudinal 
space box S — 50rs x 5rs with the number of Kapral cells 
Ncei = 1200 x 120. This is reached by adding new particles 
inside the left space box during the diffusive spreading. At 
such a density pr| — 60/25 and such a value Ncei^ the local 
diffusion rate at Q./2uj ~ 1.1 is D/Dq ~ 0.01 and during the 
time Q.st — 6.28 x 10* the diffusion propagates on a distance 
X ~ V Dt ~ 25rs that is in a good agreement with the data 
of Fig. [S] 

The situation is drastically different in presence of the 
frequency gradient g = 0.002 when the particles enter inside 
the synchronization window ^/2lu — 1 < 0.05 as it is shown 
in Fig. 15] Even if the total computation time here is 10 times 
larger than in Fig.[5]the front propagation becomes very slow 
around x/rs ~ 17 since diffusion drops strongly inside the 
synchronization window going down to very small by finite 
value D/Do ~ 3 x 10~® inside the left space box 5rs x 5rs 
with 60 particles. 

Of course, the value of the gradient chosen here is 
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Figure 7. Dependence of density of particles (arbitrary units in 
logarithmic scale) on position in the ring x/rg for e = 0.6 at the 
final time Qgt = 6.28 X 10* (red/gray for parameters of Fig. [Sjl 
and fist = 6.28 X 10^ (blue/black for parameters of Fig. |6]l. 

much larger than its real value g ~ 3rs/2xs ~ 15m/1.17 x 
lO^km ~ 10~^. Such small values of the gradient are not 
accessible for nowadays computer simulations. However, a 
more smooth, adiabatic variation of the orbital frequency 
on a scale of an epicycle radius should make the picture of 
the diffusive shock wave moving inside the synchronization 
domain to be even better justified. 

The density p{x) dependence on i at a final moment of 
time (averaged over 1% of total time) is shown in Fig. [7] for 
the cases of Figs. 15161 In the case of synchronization (Fig. [6} 
there is a sharp drop of density on a scale of 2 epicycle ra- 
dius rs- In the non-synchronized case (Fig. [5]) the decrease 
of density goes in a more smooth way (on a scale of about 
lOrs). The difference of scales in two cases is not so large 
since in both cases the diffusion is zero in the region with- 
out particles. However, the most important difference is that 
the rapid diffusive propagation goes unlimitedly in the non- 
synchronized case while inside the synchronization window 
the propagation front moves very slowly with the formation 
of sharp density drop on a scale of about 2rs- The front stays 
the longest time at the place where the diffusion is minimal 
that corresponds to the center of the synchronization win- 
dow at n = 2lu. For the value of rs oc 1/vep « 36m given 
above this gives the size of the edge Ax^ ~ 2rs ~ 70m. This 
value is in a good agreement with the observation data which 
give Axe ~ 10m especially if we take into account that the 
average value Vep is known only by an order of magnitude. 



4 DISCUSSION 

Our studies based on the SYNC model of dynamics in the 
rings of Saturn show the emergence of synchronization in 
the vicinity of the outer B ring edge. Like the Maxwell de- 
mon this synchronization makes the epicycle motion of par- 
ticles to be synchronous that practically eliminates collisions 
between them. This gives a suppression of diffusion by or- 
ders of magnitude and a formation of a diffusive shock wave 
slowly propagating inside the synchronization domain. The 
size of this front or the edge of the ring is of the order of 
a few epicycle radius being of about of a few tens of me- 
ters. This is in agreement with the present observation data 
which give its size to be about 10 meters. The size of the 



synchronization domain created by Mimas is of the order 
of 4000A:m (remarkably, the total number of synchronized 
particles, which can be estimated as > lO^'^, is huge). The 
front moves most slowly in the center of the synchronization 
domain so that it is most probable to observe it at the exact 
resonance 2:1 between the edge frequency and the frequency 
of Mimas. The observations show that the actual position is 
close to this value with a relative accuracy of 10"'' - 10"^ 

Similar synchronization effects should exist at other res- 
onances with other satellites. Our preliminary data give a 
similar picture for the outer A ring edge which is closely 
located to the 7:6 resonance with Janus (here the dimen- 
sionless kick strength is e « 1.19 but the relative frequency 
size of the synchronization window is approximately twice 
smaller due to higher order of the resonance). 

The synchronization mechanism proposed here can also 
be responsible for existence of very narrow planetary rings. 
Indeed, if initially particles are distributed inside the reso- 
nance then those which are inside the synchronization win- 
dow will remain there practically forever since the collision 
induced diffusion is switched off, while those outside of the 
window will diffuse away living a narrow ring of particles 
inside the synchronization window. 

The SYNC model used in this paper is based on several 
significant simplifications of the real particle dynamics. The 
use of the mesoscopic Kapral method for incorporating the 
collisions appears to imply no additional physical mecha- 
nisms, but the basic physics of the collisions remains hidden 
inside the parameters of the Kapral method - the number 
of the cells and the frequency of the reshuffling of velocities. 
To make the calculations more realistic, one needs to incor- 
porate realistic collision models with reasonable mechanical 
properties of the particles, to include a distribution of the 
sizes, etc. Such an extension goes far beyond the scope of 
our preliminary study. Thus we can hardly make direct pre- 
dictions, in particular, compare properties of A and B ring 
edges. 

Moreover, it appears that another simplification of our 
computational model - the use of the Nose-Hoover thermo- 
stat - is less "repairable" and requires a further justification. 
Indeed, if there is some average dissipation due to particle 
collisions, this justifies the use of the thermostat for model- 
ing the saturation of the energy pumped to the system due 
to shear. On the other hand, in our setup in the synchro- 
nized state collisions are rare so at the first glance there is 
no mechanism for energy saturation. Here it is important 
to mention, that we restricted our analysis to an ensem- 
ble of identical particles. For a particular situation in the 
outer B ring this means that we consider only large parti- 
cles that have low characteristic epicyclic velocities. If there 
is a whole range of particles of different sizes, they are all 
subject to resonant kicks by Mimas, but because their char- 
acteristic epicyclic velocities are different, the effect of the 
kicks is also different. If we assume rough equipartition of 
epicylic kinetic energies, then smaller particles have larger 
velocities, therefore for them the effective forcing parameter 
e is smaller, and as a result they only weakly synchronize or 
not synchronize (if they lie near the bottom of Fig. [2]). Col- 
lisions with these randomly moving particles may provide 
an additional dissipation that justifies the use of the Nose- 
Hoover thermostat. Further investigations (which however 
go far beyond the scope of this paper) of different dissi- 
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pation effects influencing the energy balance are needed to 
clarify this issue. 

Another simplification made - a consideration of a rel- 
atively small box of particles - is also crucial. Indeed, parts 
of the ring at different angular coordinates are statistically 
equivalent, but the dynamical equivalence is broken by the 
influence of Mimas, as the latter kicks the particles at dif- 
ferent phases. Thus the synchrony of the velocities may be 
only local, i.e. the velocities of particles are synchronized 
with the local phase of the Mimas, but this phase changes 
gradually along the ring angular coordinate. Because of the 
differential rotation, the particles synchronized at different 
phases will mix, but this effect is not taken into account in 
the box model, it has to be addressed in future studies. 

In this respect the important issue is that of how the 
synchronization of particles' velocities could be tested ex- 
perimentally. Indeed, the theory above predicts a drastic 
narrowing of the width of the velocities distribution, cf. 
Fig. [1] However, because the velocities themselves are very 
small, there is no much hope to observe the distribution 
of them directly, e.g. by Doppler measurements or image 
analysis. Thus one has to rely on implicit observations. For 
example, one may expect that the adjustment of veloci- 
ties influences other dynamical feat ures recently observed i n 
the rings, like propeller structures ijSremcevic et al.l l|2007t l: 
IS pahn fc Schmidd ()2006l l ) . However, for such an implicit test 
one has to find such structures quite near to the sharp edges 
described above. Remarkably, near the outer edge of ring 
B a scrambled pattern is observed in recent Cassini images 
( see URL: http://saturn.jpl.nasa.gov/photos/imagedetails/ 
index. cfm?imageld=2984 ) that is probably due to a 
gravitational clumping of particle (cf. simulations by 
iLewis fc Steward (|2005l ) of a similar structure at the Encke 
gap). We can speculate that the synchronization of large 
particles would generally enhance their ability to clump due 
to gravitational and adhesive forces, because their relative 
collision velocities become rather small. On the other hand, 
if smaller, non-synchronized particles, leave the ring and en- 
ter the gap, where large particles are missing, their diffusion 
will be reduced as well, because the collisions between small 
particles are rare. Observation of such a gradient in the size 
distribution of particles near the edge may also be inter- 
preted as an implicit confirmation of the suggested mecha- 
nism. 

We note that in the literature several mechanisms for 
explanation of sharp edges of rings have been discussed. One 
process to form and maintain a sharpness is a shear rever- 
sal induced by a strong gravitational perturbation near a 
resonance with a moon. In the unperturbed ring Keplerian 
shear leads to an outward viscous transport of angular mo- 
mentum. The distortions induced by a perturber may locally 
reverse the shear and, when integrating over the azimuth, to 
a net an gular momentum transport that i s directed radiall y 
inward (jBorderis et all Il982l . 1 1981 1 19891 : iMosaueiral 1 19961 ) . 
This process can lead to an inward migration of material 
from the close vicinity of the resonance location with the 
moon, and in this way clear a gap at this radial position in 
the ring. Recently Lewis et al. reported on a "negative dif- 
fusion" mech anism, where the parti cles migrate to areas of 
high density l|Leezer fc Lewis! (|2006l )). In their calculations, 
however, only a single pass by the moon was simulated so 
no resonant effects, like those we focus on in this study, were 



possible. Probably, further observations, together with fur- 
ther extensive numerical simulations, will allow one to test 
which mechanism is responsible for the sharp edges of the 
rings. It appears highly desirable to compare various setups 
in numerical simulations as well. 

Finally, we note that, as discussed by 
IChepelianskii et al.l ()2007h . the elimination of colli- 
sional diffusion may appear also for char ged partic l es in a 
magnetic fie l d like 2D electron gas (see iMani et al.l l|2002h ; 
IZudov et al.l ll2003h) and electron and ion clouds in a Paul 
trap (see, e.g., Mortenson et al.l (|2006l )). Due to that it can 
be rather interestin g to try to m a ke lab oratory experiments 
with traps (see e.g. iMaior et al.l |2005f )) which would allow 
to model rings of Saturn in laboratory experiments with 
cold ions. 
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diffusion constant for P, because in the absence of collisions 
this quantity is exactly conserved. 



APPENDIX B: DERIVATION OF MIMAS'S 
KICK FORCE 

Here we derive the strength of kick force produced by Mi- 
mas on the epicycle motion of particles inside the Saturn 
ring B. Consider the effect of a gravitational action of the 
Mimas having mass mm and semi-major axis am on a parti- 
cle having axis a (both orbits are nearly circular). In a frame 
fixed with the particle, the graviational acceleration in the 
outer direction is Vx = GMmd~^ cosi/; where d and ip is the 
distance and the angle from the particle to the Mimas. De- 
noting the angle from Saturn to Mimas as (/> = (f2 — ^m)t 
where Q.m,^ are Kepler frequencies, we can write cos?/) = 
d~ ""^(cim cos — a) and d^ = (omsin^)^ -I- {amcoscj) — i)^- 
The total change of the velocity component due to the 
force by Mimas is given by the integral J dotVx ■ We attribute 
this change to a "kick", because the force is non-zero only 
when Mimas is close to the particle. The resulting expression 
is 

^ Gmm. r r^(cos0-r) dcj) 

(J7-n™)a2 (l + r2- 2r cos 0)3/2 ^ > 

This calculation gives Avx ~ 3.2 ■ 10~^m/s. When normal- 
izing by a typical epicycle velocity 5 • 10~^m/s we get a 
numerical value of our parameter e « 0.64. 



APPENDIX A: THE HAMILTONIAN FORM OF 
THE HILL EQUATIONS 

Here we descri be the Hamiltonian form of the Hill equations. 
According to iStewarti (199 1), the Hill equations ^ can be 
viewed as a Hamiltonian system with 

H ^ ^[{px + nyf + {Py-Qxf]-^Q^x''-xF4t)/mp (Al) 

For the simul ation and the a nalysis another representation, 
also given bvlSt ewartI (jl99ll ). where the epicyclic and shear 
motion are effectively separated, is more convenient. One 
introduces canonical variables /, (/>, P, Q according to 





In these variables the Hamilton function reads 
3 ^2 / 2 „ fW 



H = ni-^P'-\^-P-^ -cos <j,]Fx{t)/m, (A4) 

Canonically conjugated variables /, (j) are the action-angle 
variables for the epicyclic motion. When introducing the 
Nose- Hoover thermostat, we adjust variable I only, mod- 
elling in this way the balance of this part of the total energy. 
Canonically conjugated variables P, Q describe the shear, 
the conserved quantity P corresponds to the conservation 
of the angular momentum. Note that according to (|A2|) . 
variables P and Q can be viewed as "center of mass" coor- 
dinates for the rotating particle, it is especially convenient 
to calculate the diffusion rate in a;-direction in terms of the 
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